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ABSTRACT 


A  method  is  developed  for  calculating  natural 
frequencies  of  cracked  rectangular  plates  and  of  cracked 
cylindrical  shells  of  rectangular  planform.  In  either  case 
the  crack  is  rectilinear  of  arbitrary  length  and  of  arbitrary 
location.  The  analysis  is  based  on  finite  Fourier  transform 
method  of  discontinuous  functions  in  the  form  suggested 
earlier  by  the  author  (Refs.  1,2,3,6,24).  For  the  shell  the 
problem  is  described  by  Donnell's  equations.  For  either  case 
the  discontinuities  of  the  displacement  and  of  the  slope 
across  the  crack  are  the  unknowns  of  the  problem.  The 
unknown  amplitudes  of  those  discontinuities  are  determined  by 
satis-  fying  the  boundary  conditions  at  the  crack's  edge  (2 
for  the  plate  and  4  for  the  shell) .  This  requires 
differentiation  of  Fourier  series  representing  discontinuous 
functions  and  is  achieved  by  using  generalized  Green-Gauss 
theorem.  The  square-root  singularities  of  the  bending  moment 
(and  also  of  the  normal  and  shear  stresses  for  the  shell)  are 
built  into  the  solution.  Ultimately  the  problem  is  reduced 
to  an  infinite  system  of  linear  algebraic  equations.  The 
method  of  reduction  is  applied  to  the  characteristic 
determinant  of  the  problem  and  numerical  values  of  the  lowest 
two  frequencies  are  found,  as  examples,  for  (1)  a  simply 
supported  square  plate  with  symmetrically  located  crack, 
parallel  to  one  edge  for  various  d/a  ratios  (d-length  of  the 
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crack,  a-length  of  the  plate  side);  (2)  for  a  simply 
supported  square  plate  with  a  crack  of  fixed  length,  parallel 
to  one  edge  at  varying  distances  from  it;  (3)  for  a  simply 
supported  square  plate  with  a  diagonal  crack  of  varying 
length;  (4)  for  a  simply  supported  cylindrical  shell  of  a 
square  planform  with  a  crack  of  varying  length  located  at  the 
apex  parallel  to  the  straight  boundary  for  various 
radius/edge  length  ratios. 
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For  the  past  two  years,  the  author  has  carried  out  a 
research  program  for  the  Department  of  the  Navy  that  was 
aimed  at  the  analysis  of  natural  vibration  of  cracked 
rectangular  plates  and  of  cracked  shallow  cylindrical  shells 
of  rectangular  planform.  The  results  of  this  investigation 
have  been  reported  elsewhere, [1] ,  [2],  [3],  here  more  details 
of  derivation  are  added. 

The  report  consists  basically  of  two  parts:  in  Part  1 
the  cracked  plate  is  analyzed  while  the  cracked  shell  is 
investigated  in  Part  2. 

The  studies  of  bending  and  vibration  of  finite  cracked 
plates  are  limited  to  few  papers.  Bending  of  a  cracked 
rectangular  plate  was  first  investigated  by  Keer  and  Sve  in 
[4].  Their  analysis  was  limited  to  such  location  of  the 
crack  that  allowed  reduction  of  the  problem  to  a  dual  series 
equation.  Hence  the  crack  was  confined  to  a  position  along 
the  symmetry  axis.  Analogous  method  was  applied  by  Stahl  and 
Keer  [5]  for  analysis  of  natural  vibration  and  stability  of 
rectangular  plate  and  was  bounded  by  similar  limitations  as 
encountered  in  [4].  In  [6]  Solecki  attempted  to  remove 
existing  restrictions  by  developing  a  method  that  would  allow 
to  study  rectangular  plates  with  arbitrarily  located  crack. 

He  partially  suceeded  by  developing  a  method  based  on  the 
combination  of  finite  Fourier  transformation  and  of  the 
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generalized  Green-Gauss  theorem.  A  cracked  rectangular  plate 
was  discussed  as  one  of  the  examples.  Numerical  data  were 
not  however  obtained  partially  because  the  singularity  of  the 
curvature  at  the  tips  was  not  explicitly  isolated.  The  only 
other  paper  dealing  with  a  similar  topic  was  written  by  Ali 
and  Atwani  [7]  who  used  a  version  of  Rayleigh's  method  to 
study  natural  vibrations  of  rectangular  plates  with  cutouts. 

Numerous  papers  were  devoted  to  analyses  of  vibrations 
of  thin  shells  of  various  geometries  (see,  for  example,  Dym 
[8]  and  Kraus  [9]  for  references) .  All  reported  studies 
however  were  applied  to  homogeneous  (uncracked)  shells. 
Apparently  no  investigation  of  vibration  of  cracked  shells 
was  done  so  far.  Considerable  mathematical  difficulties  are 
encountered  when  attempting  to  solve  this  problem  formulated 
by  a  coupled  system  of  partial  differential  equations 
describing  functions  with  discontinuous  derivatives. 

Let  us  survey  briefly,  related  studies  on  statics  of 
cracked  cylindrical  shells  noting  that  almost  all  of  the 
analyses  refer  to  infinite  shells.  The  first  attempt  to 
determine  deflection,  stresses  and  stress  intensity  factor  of 
a  cracked  shallow  cylindrical  shell  is  due  to  Folias  [10]  who 
reduced  Marguerre-Reissner  differential  equations  describing 
the  problem  to  a  system  of  dual  integral  equations.  This  in 
turn  can  be  transformed  to  a  set  of  coupled  singular  integral 
equations  of  the  Cauchy  type. 


This  technique  of  reducing  a  system  of  partial  dif¬ 
ferential  equations  to  a  system  of  dual  integral  equations, 
and  subsequently  to  a  system  of  singular  equations,  seer.j  to 
be  the  favorite  of  most  authors  of  papers  on  statics  of 
cracked  shells. 

A  similar  type  of  shell  to  that  discussed  in  [10] ,  was 
treated  with  similar  mathematical  devices  by  Copley  and 
Sanders  in  [11].  Erdogan  and  Kibler  [12]  also  studied 
shallow  cylindrical  shells  using  Marguerre-Reissner  equations 
and  reducing  them  to  a  system  of  singular  integral  equations. 
Their  solutions  however  is  valid,  unlike  Folias's  solution 
[10],  also  for  relatively  long  cracks.  Keer  and  Watts  [13] 
based  their  analysis  of  a  complete  cylindrical  shell  with  a 
circumferential  crack,  on  the  equations  of  three-dimensional 
theory  of  elasticity.  Sanders  [14]  used  his  equations  to 
analyze  an  infinite  cylindrical  shell  with  a  circumferential 
through-crack. 

Few  solutions  exist  that  are  based  on  10th  order  shallow 
shell  theories  (Naghdi  equations) .  These  are  also  obtained 
by  reducing  the  differential  equations  to  singular  integral 
equations  via  dual  integral  equations.  Krenk  [15]  used  this 
theory  (which  takes  into  account  transverse  shear 
deformation)  to  analyze  cylindrical  shell  with  an  axial 
crack.  Most  recently,  Delale  and  Erdogan  [16]  also 
investigated  the  effect  of  transverse  shear  on  the  behavior 
of  circumferentially  cracked  cylindrical  shell. 
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In  [17]  Delale  analyzed,  among  others,  also  an  axially 
or  circumferentially  cracked  infinite  cylindrical  shell.  An 
original  approach  often  applicable  to  quite  general 
geometries,  partially  based  on  Sanders  paper  [18] , 
characterizes  works  of  Simmonds  and  colleagues  (see,  for 
instance,  [19],  [20]  and  [21]}.  Finally,  still  another 
technique  due  to  Lakshminararyana  and  Murthy  [22]  and  [23]  is 
valid  for  short  arbitarily  located  cracks  in  infinite  shallow 
shells,  it  is  based  on  formulating  the  problem  in  elliptic 
coordinates.  This  allows  the  crack  edges  to  become 
coordinate  lines  (as  a  limit  of  an  elliptic  cutout) .  This 
procedure  although  elegant  and  useful  in  establishing  close 
form  expressions  for  local  quantities  would  be  too  cumber¬ 
some,  even  if  applicable,  to  deal  with  global  quantities, 
because  it  would  require  manipulating  infinite  series  of 
Nathieu  functions. 

In  the  static  case  differential  equations  describing 
various  shell  theories  are  amenable,  for  certain  geometries 
and  special  crack's  locations,  to  a  system  of  dual  integral 
equations.  These  are  in  turn  transformed  into  singular 
integral  equations  which  are  usually  solved  numerically. 

When  inertia  terms  are  included  in  the  differential  equations 
and  when  in  addition  the  crack's  location  and  its  relative 
length  are  arbitrary,  the  above  procedure  encounters 
considerable  difficulties  or  is  even  not  applicable  (in 
particular  for  finite  geometries) . 


In  part  1  of  the  present  report  a  method  is  demonstrated 
for  determining  frequencies  of  steady-state  vibration  of  a 
rectanqgular ,  simply  supported,  isotropic  plate  with  an 
arbitrarily  located  crack  (Fig.  1) .  It  is  based  on 
application  of  finite  Fourier  sine  transformation  in 
conjunction  with  the  generalized  Gre^n-Gauss  Theorem.  This 
method,  as  was  shown  in  [6]  and  [24],  eliminates  considerable 
amount  of  tedious  integration  by  parts. 

After  applying  double  finite  Fourier  transformation  to 
the  differential  equation  governing  the  problem,  and  after 
using  the  inversion  theorem  one  obtains,  as  usual,  a  system 
of  integral  equations  with  respect  to  the  unknown 
discontinuties  of  the  deflection  and  of  the  slope  across  the 
crack.  It  is  known  [25]  ,  [1]  that  using  as  the  unknowns 
higher  order  derivatives  one  improves  the  convergence  of  the 
resulting  infinite  series.  Therefore  using  the  condition 
that  the  bending  moment  at  the  crack  equals  zero  one  can 
replace  the  integral  involving  unknown  discontinuity  of  the 
deflection  by  the  discontinuity  of  the  curvature  in  the 
direction  normal  to  the  crack.  The  unknown  quantities  are  in 
turn  expanded  into  Fourier  series.  Since  the  curvature  is 
square-root  singular  at  the  tips  therefore  it  is  represented 
as  a  sum  of  a  regular  part  (expanded  into  Fourier  series)  and 
of  a  singular  part  with  known  strengths  of  the  singularities 
multiplied  by  unknown  coefficients  determined  from  the 
conditions  of  the  problem.  Finally,  the  unknown  infinite 
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sequences  of  Fourier  coefficients  are  obtained  by  applying 
the  conditions  that  the  bending  moment  and  the  shear  force 
across  the  crack  are  zero.  Application  of  these  conditions 
requires  dif f erentation  of  Fourier  series  of  discontinuous 
functions.  A  novel  approach  is  suggested  in  relation  to 
this  operation.  Determining  Fourier  coefficients  of 
derivatives  constitutes  in  this  case  a  major  undertaking. 

This  effor  is  reduced  to  a  fraction  by  applying  again  the 
generalized  Green-Gauss  theorem.  In  Part  2  of  the  report 
Donnell  equations  are  used  to  describe  vibration  of 
cylindrical  shells.  The  same  method  as  in  Part  1  is  applied 
again  here.  Its  main  feature  consists  in  expanding 
displacement  functions  discontinuous  inside  the  region  of 
interest  into  infinite  series  and  then  proceeding  according 
to  the  rules  of  integral  transformations  developed  previously 
for  discontinuous  functions  (application  of  integration  by 
parts,  Green-Gauss  theorem,  etc.).  Differential  equations 
are  transformed  as  the  result  of  these  operations,  into 
algebraic  equations  with  respect  to  the  transforms  of  the 
unknown  displacements.  Application  of  the  inversion  theorem 
allows  to  determine  formally  the  displacements  that  still 
include  at  this  stage  certain  integrals  depending  on  unknown 
quantities.  Specifically,  the  unknown  functions  appearing  in 
the  integrands  are  the  quantities  discontinuous  across  the 
crack  (discontinuities  of  the  first  kind) . 


In  addition,  depending  on  the  boundary  conditions  on  the 
shell's  circumference,  other  unknown  integrals  may  appear. 
These  include  certain  unknown  boundary  quantities.  The 
objective  of  the  next  step  is  evaluation  of  the  unknown 
functions  in  the  integrands. 

At  this  point  it  is  very  important  to  realize  (1)  which 
unknown  functions  are  singular  at  the  ends  of  their  range  of 
definition,  and  (2)  what  are  the  strengths  of  singularities. 

Sufficient  experience  and  information  was  gathered  by 
various  researchers  to  answer  these  questions.  It  is 
therefore  possible  to  represent  the  unknown  functions  as  the 
sum  of  a  regular  part  in  the  form  of  infinite  series  with 
unknown  Fourier  coefficients,  and  of  a  singular  part  of 
known  form  multiplied  by  an  unknown  constant  (the  intensity 
factor) . 

Subsequently  available  boundary  conditions  at  the  edge 
of  the  crack  are  applied.  At  this  stage  particular  care  must 
be  exercised  when  differentiating  infinite  series 
representing  functions  with  "inside"  discontinuities:  they 
never  can  be  differentiated  term  by  term.  The  appropriate 
formulae  for  differentiation  are  derived  applying  again 
generalized  Green-Gauss  theorem.  Once  all  these  mathematical 
obstacles  are  overcome  the  final  result  is  obtained  in  the 
form  of  an  infinite  system  of  linear  algebraic  equations.  In 
the  case  of  natural  vibration  the  system  is  homogeneous  and 
depending  on  the  unknown  natural  frequencies.  Thus  the 


characteristic  equation  is  in  the  form  of  an  infinite 
determinant  equated  to  zero.  As  usual,  assuming  that  the 
system  is  regular,  one  applies  the  method  of  reduction 
obtaining  frequencies  with  any  desired  accuracy.  Adaptation 
of  this  procedure  to  higher  order  theory  would  not  present 


any  major  difficulties. 


II. 


Vibration  of  Cracked  Rectangular  Plates 


1.  Series  Solution  of  the  Differential  Equation 


The  differential  equation  which  governs  the  amplitude  of 
flexural  vibration  w(x,y)  of  the  plate  shown  in  Fig.  1  is 

D  V2V2w-pw2w  =  q(x,y)  (1) 

with  the  boundary  conditions  of  simple  support  around  the 
contour  and  the  conditions  of  vanishing  bending  moment  Mu  and 
vanishing  shear  force  Qu  at  the  crack's  edge  where  u  is  the 
direction  normal  to  the  crack. 

The  solution  of  Eq.  (1)  is  assumed  in  the  form 

SlnV  Sin6ny  <2> 

In  order  of  finding  wmn  one  applies  finite  Fourier  sine 
transformation  to  eq.  (1).  As  in  [24]  this  is  performed  in 
conjunction  with  the  generalized  Green-Gauss  theorem 

j|  F,.  dS  =  |  Fn.  dr  +  |  Fn.  df  (1=1 ,2,3)  (3) 

S  f 

where  T  is  the  outer  contour  of  the  plate,  F  -  its  inner 
contour,  n^  and  n^  are  the  components  of  unit  vectors  normal 
to  T  and  F  respectively  and  where  F  is  any  tensor. 


Application  of  this  procedure  and  subsequent  solution  of 


Thorne  [25]  is  applied,  that  is  the  unknown  Aw  is  replaced  by 
the  unknown  discontinuity  of  the  curvature  in  the  direction 
normal  to  the  crack:  A  32w/3u2.  This  is  achieved  by 
utilizing  the  condition  that 


0 


(8) 


Integrating  the  first  integral  in  (4)  by  parts  and  using  (8) 
one  obtains  the  following  relations 


o 


(9) 

The  unknown  slope  discontinuity  is  now  expanded  into  Fourier 
cosine  series 

A  lu  *  f  2  W  cos  IT  »  ce(0,d)  (10) 

m 

The  unknown  discontinuity  of  the  curvature  has  square-root 
singularities  at  the  tips  of  the  crack.  It  is  therefore 


represented  as  a  sum  of  a  the  unknown  regular  function  which 
is  expanded  into  Fourier  sine  series  and  of  a  singular 
function  with  square-root  singularities  at  the  end  of  the 


interval  ce  (0,d)  and  adjusted  in  such  a  way  that,  in  concert 


with  the  physical  conditions  of  the  problem,  both  Aw  and 


A3w/3u  vanish  at  the  ends. 


This  leads  to  the  following  representation: 


2  “  u  e,  kirc  A  3  “  Wk  .  1  .  1  » 

-r  E  W.  sin  —r~  +  -  E  T—  ( - +  - ) 

d  k=1  k  d  v  VS  k=l  k  /E 


—  E  ^  - 

_  nr  i-_i  ->  K  rz 


n/S  kxl  ,3  K  SE  /d-c 


Presently  eqs.  (10)  and  (11)  are  used  in  (4)  which  is  in  turn 


substituted  into  eq.  (2)  leading  to  a  formal  solution  of  the 


problem  which  still  includes  unknown  Fourier  coefficients  Uk 
and  W^. 


2.  Boundary  Conditions  and  the  Characteristic  Equations 


In  order  to  determine  two  infinite  sets  of  unknown 


Fourier  coefficients,  two  boundary  conditions  available  at 


the  crack's  edge  are  applied 


T^j-r  +  v  T~r  *  0  at  x=x^  +  c  cosif>,  y=y1  +  c  simj> 


i!"  +  (2-v)  3\ 
SiF  '  v)  3U3C3 


*  0  at  x=x-|  +  c  cosip,  y=y-j  +  c  sim/> 


•^y:v:v:£vS::v£%^ 


It  is  seen  that  the  series  (2)  must  be  differentiated  ar.d 


that  this  can  not  be  done  term-by-term  because  (2)  represents 
a  function  with  discontinuous  second  and  third  derivatives. 

As  indicated  in  [1]  proceeding  in  traditional  way  will  result 
in  proper  formulae  being  devised  but  at  a  cost  of 
considerable  amount  of  cumbersome  manipulations.  The 
procedure  devised  in  [1] ,  based  on  Green-Gauss  theorem  is 
generalized  here  to  include  the  use  of  straight  line 
discontinuity  arbitrarily  located  within  the  plate's  area. 

Let 

♦(Iran)  *  si"V  sinV 
♦(2mn)  -  C0SV  s1nV 
*(3mn)  *  si"V  C0SV 

♦(-tom)  -  C0SV  C0SV  (14) 


Assume  that  we  want  to  obtain  the  expression  for  3w/9u  where, 
as  before,  u  is  the  direction  normal  to  the  crack. 

Let  F  =  w  4>(2mn)  be  substituted  into  (3)  for  i  =  1. 

The  following  relation  results 


|J  fe  (w  *ta«n)>dS  1  |  w  *(2mn)uxllF 
S  r 


(15) 


where  the  integral  along  r  vanished  because  of  the  assumed 
boundary  conditions. 


18  - 


Expanding  the  integral  on  the  left  hand  side  of  eq.  (15) 
and  taking  into  account  that 


\\  k*(2mn)  *  ♦(lmn) 


and  that 


JJ  w  +(lmn) 


ds  *  w. 


one  obtains  finally 


f{  !r  “SV  sinV  *  Vm,  *  \  w  ♦(2n.n)“xdf 


One  can  easily  find  now  the  formula  for  3w/3x : 


I?  *  h  j,  |w(sinV!r  Vf  5,nV 


1^  l  I  (Vmn  ♦  |  w(cosV  s1n6ny)f  u>!dp  C0SV 


s1nV 


In  a  similar  way,  starting  with  the  expression  F  =  v  <b(2mn) 


one  obtains  from 


II  3y  (w  *(3mn)>  dS  *  f  w  d(3mn)  “y  df 


‘  *  -  -  •  VT4jf^  - 


the  following  formula 


!r  iFraJ,  \  w(sinv>f  uy  df  s1nv 


",  2  (6„  *  j  *<si"V  eosV).  “y  dr  sinV  C0SV 

npl  nsi  _  i 

r  (19) 


Now  using  known  expression  for  directional  derivative 


3w 

3u 


3w 

3x 


cost  +  ~ 


simp 


(20) 


yields  the  desired  result. 

The  same  procedure  was  used  to  obtain  higher  order 
derivatives  appearing  in  the  boundary  conditions  (12)  and 
(13).  Resulting  system  of  equations  still  depends  on  the 
local  coordinate  c  (along  the  crack) .  This  dependence  is 
removed  by  multiplying  the  equations  obtained  from  (12)  and 
(13)  by  cos  (inc/d)  or  sin  (Jlirc/d)  respectively  (where  1  - 
0,1,... «)  and  integrating  the  results  with  respect  to 
c  e(0,d)  . 

This  yields  an  infinite  system  of  linear  algebraic  equations 
of  the  following  form 


tM  “k  C|kt  +  k4,3  "k  ^  *  4o  “k  Cjkk  '  0>  l-“-'  - 


Wk  C4kt  + 


«  00 
?  -  c«;n>  +  ^  U.  cc, .  =  0, 


1-1.. 


where  U  =  -2  I  U.  and  where 
0  k=i  k 


Clkt  ■  v  (*a5s  .  ■  T4ffiOk  CwZ2ni0k  *  JF  c3m0- 

m- 1 


n] 


+  E  n  T,n_0  [ttZ-, 


3^r 


‘5c  "  '40nJt  lT,L]Onk  "  ~2lT  E10n^} 


ao  ao 

*h  *  S  -r-  «*,_  T^„,  ♦  X,„  T.  .H-irZ. 


V  mr1  "Almn  3mn£  A2mn  ’toi". 


3^tT 


^lmn^  ^  ^  "imn  "^Amn 0  )( 


2k  Imn7  v  3mn  3mn£  4mn  4fran£/v  2mnk 


E^.)) 


k  3mn' 


(22) 


..  C^2ti 


2kJL  ~  kv  *"*a5s  _Z,  m  T4mQ£  E4mO  "  a5r  E_  n  T/mn«  E 


m 


=1  4m0£  4m0  "  a5c  ^  n  ‘40n£  ^On 


00  oo 


+  mfi  nf]  Amn  ^Xliro»  E2mn  +  X3mn  E4mJ  T3mn£ 


^X2mn  E2mn  +  X4mn  E4mn  ^  ^4mn£^ 


(23) 


Xk  “  “ 

c3k£  =  {all  *  T4mOk  T4mo£  +  all  \  T40nk  T40n£ 
m-i  n- i 


+  2. 


m= 


1  n=1  Amn  L'  ®mn  ®mn  3mnk 


+  ^X7mn  T3mn£  +  X8mn  T4mn£^  T4mnk*  +  1T5anfi 


11°0£ 
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where  6g^  is  Kronecker's  delta,  and  =  1  -  1  fi 
The  coefficients  c^^,  c5kl#  c6kl  are  obtained  from 
c2kl  an<*  c3kl  respectively  by  making  the  substitutions  shown 
in  Table  1.  The  expressions  for  the  terms  appearing  above 
are  given  in  the  Appendix  1.  The  characteristic  equation  is 
obtained  by  equating  to  zero  the  infinite  characteristic 
determinant  of  the  system  of  linear  algebraic  equations  (21) . 
In  practice  the  method  of  reduction  is  applied  and  the 
frequencies  are  determined  from  a  finite  characteristic 
equation. 


3.  Examples  Is  A  Plate  with  a  Parallel  Crack 


As  the  first  example  consider  a  plate  with  a  crack 
parallel  to  one  edge  (see  Ref.  [1]  and  Fig.  2) .  In  the 
present  case  eq.  21  can  be  simplified  because  the  double 
series  appearing  there  (see  Eqs.  27*24)  can  be  reduced  to 
single  series  to  yield  the  following  system: 

?  kL  (Cm  f.  -  rr  Cc_)  "  ?  lii,  ^2kn  =  ®  P  =  0,2,..'° 
c-1,3  k  4pk  P  5p  k=2,4  k  Zkp 


3  «%  -  f  V  \.f>4  V 0  p”-3-" 


The  dimensionless  quantities  appearing  in  eq.  (25)  are  given 
below: 


npk  '  Pp  m4i3  ~iF^  [U-v)V  -  f!]Zim(e.e)  p  •  1.3.. .;k 


2kp  *  Pm=f>3  (C(l-v!)m'  -  f’lZJe,') 


*  (l-v)[Z^(0.0)  -  Z,(2e,0)]> 


p  -  0,2,. ..;k  =  2,4,. 


a _  a. 


c4pk  -  aV  f  -® 

P  *  m-1 ,3  m 


C,  -jfi  f  1 
5P  aw  _ 


[(l-v)2m4  -  f2]Z.  (e,e)  p  =  o,2,..;k  - 


J.2  rfV”  [C0S  ^  +  s1n  T* 


[(l-v)2m4  -  f2]Z  (e,e) 


P  =  0,2,. 


Pkp  =  Pip  J^  3  Pmmikm  {_  0-v)m2  +  2(l-v)m2(f2  -  ni-JZ^e.e) 


+  Cf2  -  OV)«»]L  (e,e)> 


^  »P  B  1,3,... 


RP  ’  *»P  «r-?,3  iHMiF  [C0S  T4  +  s'n  ^  S(«J)]{  -  (l-v)m! 

*  ZO-'>)'"2(fJ  -  ■-)  Z5m(e,e)  *  [f!  -  (l-u2)m"]z3m(e,e)J  p 
/where  ♦  =  d/a,  y  *  a/b,  and  e  =  e^/b 

Z.  (£..£,)  =  y  ?  n  sin.mrgi  cos  nirg, 

Y  n*l  1PP  +  m2)2  -  'f*  » 

7.  (r.  r.\  =  9  sin  nir£i  sin  nwE, 

^mV4l**2/  n*l  TyV  +  ip)'  -  }* 


7.  (c.  e.\  -  ..2  9  n2  cos  n 

Y  ni}  7yV  ♦ 


nir£i  cos  mr£,  . 
+  m*)5  -  f2 


1  raft  a»?h  • 


The  characteristic  equation  is  obtained  by  equating  to  zero 
the  finite  characteristic  determinant  of  the  system  (25).  in 
practice  method  of  reduction  is  used  to  obtain  the  natural 
frequencies,  that  is  the  number  of  rows  and  columns  is 
limited  to  a  certain  finite  number. 

For  a  square  plate,  Y=l,  and  the  functions  (27)  become: 


zi  (Ci,C2)  =  hr  {— ~-n  +  sinh 

8f  sinh  tr 

-  C  S1nh  +  sinh  FTr(l-Ei-E?)itf +f  ] 

sinh  Th42+f 

where  C  =  f  1  for  >  ** 
for  Ci  <  €2 

Z~(Ci,i;2)  =  h  {  -c-°-sh  ^O-ki-CzDvffi7^-  cosh 
8f  sinh 

-  cosh  [tt(1-Ui-€2  1  +  cosh  [Tr(l-Ci-E,)»42+f1  j 

sinh  Tr/m^+f 

Z*.Ui.€a)  =  h-  cosh  [ird-k^CzD/i^  ±  cosh  [tt(1 -Ei-^).'^  ] 

*  8f  sinh 

+  rfF+f  cosh  [nd-jgi-gzl  )y|iri7tf  +  cosh  l>(l-Ei-£,)/mz+f  1  } 

sinh  ir/mz+f 

l^Uuit)  =  ^  {  *-CP.sh  Cir(1-l€i~ga I )w1ng^F  ]  +  ir  cosh  [ff(l-C,-C,)/Ff  1 

2/mz-f  sinh 

_  1  cosh  [tt(1-1€i-C2  I  )/mz+f  1  +  tt  cosh  [TrO-Ei-Esl^nF+f  1 

2/mz+f  sinh  7r/mz+f 

2f  . 

‘  ^r-fT  i  (28) 
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The  system  of  algebraic  equations  (25)  was  also  put  in  more 
compact  form  (see  [26])  by  introducing  the  following 
definitions: 

Wk  for  k=l,3,... 

Ak  *  { 

for  k=2,4,... 

C4pk  “  P* C5p  for  p=0,2,... 

^2kp  ^or  k-2,4,...;  p=0,2,... 

ckp  *  {  ] 

Pkp  _  W  Rp  for  p=l  ,3,. . . 

Clplc  for  k=2,4,. . . ;  p=l,3,... 

Now  the  system  of  eqs.  (25)  takes  the  form 

kf,\  Ckp  *  0  P=0,1,..- 

and  the  characteristic  equation  becomes 

det|Ckp|  =  0  k=l,2,..«>;  p=0,l..«  (3i) 

Characteristic  equation  (31)  was  first  solved  for  the  crack 
located  on  the  symmetry  line  y  *  b/2  =  a/2  (see  fig.  2) . 

This  allowed  for  commparison  with  the  existing  results  given 
in  [5].  Table  2  presents  the  smallest  values  of  the 
frequency  factors  f  corresponding  to  symmetric-symmetric  and 
symmetric-antisymmetric  ("opened  crack")  vibrations  and 


(29) 


(30) 
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compares  them  with  values  calculated  in  [5].  It  should  be 
noted  that  the  present  results  were  obtained  by  considering  a 
20  X  20  matrix  (in  few  cases  a  30  X  30  matrix).  For 
"opened-crack"  vibrations  the  results  were  not  sufficiently 
accurate  for  a  crack  extending  almost  to  the  edges. 

When  crack  is  not  located  on  the  symmetry  line  both 
modes  are  coupled.  The  frequency  factors  f  for  2  lowest 
modes  were  calculated  for  d/a  B  0.5  and  for  various  relative 
distances  of  the  crack  from  the  edge  (0  <  e^/a  <  0.5). 

Results  are  given  in  table  3  and  represented  graphically  in 
Fig.  3.  It  is  interesting  to  note  that  the  smallest  value  of 
the  frequency  factor  corresponding  to  the  second  mode  is 
obtained  for  a  crack  located  approximately  at  e  /a^O. 25. 


Example  2:  A  Diagonally  Cracked  Square  Plate 


Ckp  =  C1 (2k-l )(n/2)  +  c2(2k-l ) ( p/2) *  k=1»5*9.---;  P=0,2,... 

Ckp  =  C1 (2k-l ) (p/2)  ’  k=3*7.H p=0,2, . . . 

Ckp  ~  c3(2k) (p/2)  *  k=2«3»-**;  P=0,2,... 

Ckp  =  c4(2k-1)(2p-l)/2  +  c5(2k-l )(2p-l )/2  *  k=l  ,5,9.. ;  p=l,3,... 

Ckp  =  c4(2k-l)(2o-l)/2  +  k=3’7*n »••• J  Psl,3,... 

Ckp  =  c6(2k)(2p-l )/2  »  ks2»4.---;  P=l,3 

Characteristic  equation  for  the  system  (32)  becomes 

det  |  |  =  0  k=l,2,...«;  n=0,l,...» 

The  geometry  of  the  plate  and  of  the  crack  allows  for  modes 
which  are  symmetric  or  antisymmetric  with  respect  to  the 
diagonal. 

For  such  modes  the  natural  frequencies  have  been  found 
for  various  relative  lengths  of  the  crack:  C  =  d/a.  The 
results  are  presented  in  Fig.  4.  When  the  crack  extends  to 
the  corners  the  corresponding  frequencies  should  equal  these 
of  a  right  triangular,  equilateral  plate  simply  supported 
along  the  legs  and  with  the  hypothenuse  either  free 
(symmetric  vibration  of  the  cracked  plate)  or  "sliding" 
(antisymmetric  vibration  of  the  cracked  plate) . 

Unfortunately  these  data  are  not  available  for  this  type  of 
supports  of  a  triangular  plate.  The  indirect  verification  of 
the  results  is  obtained  by  noting  that  the  frequencies 


of  both  modes  seem  to  converge  to  the  same  value  when  C  is 
increasing  to  ^2.  It  was  not  possible  to  determine  the  exact 
value  of  this  limit  because  the  convergence  of  the  series 
involved  decreases  when  C  approaches  /if.  This  pheonomenon  is 
similar  to  the  one  noted  in  ref.  [24]. 


Ill 


Vibration  of  Cracked  Cylindrical  Shells 


1.  Series  solution  of  the  differential  equations 


Consider  a  cracked  shell  shown  in  Fig.  5.  It  is  assumed 
for  simplification  that  the  deflections,  normal  displace¬ 
ments,  bending  moments  and  shear  stress  vanish  on  the 
contour.  Validity  of  Donnell  Equations  [8]  is  assumed: 


(i4 

S> 

32v 

3C3n 

♦ « Hf  ■  -<//c 

(1) 

l+v 

~r 

32u 

3£3n 

<¥ 

14  + 
W 

£*  ♦  fc  ■  V/e 

(2) 

v k 

♦  !* 
3n 

+  a 

+  k2  V2V2)w  -  -qnR2/C 

(3) 

where  v  is  Poisson's  ration,  k  is  the  thickness  parameter, 
and  C  is  the  extensional  rigidity 

k2  -  (£)2  .  C  *  Eh/(l-v2)  (4) 

where  R  is  the  radius  of  curvature,  h  is  the  thickness  of 
the  shell  and  where  the  body  forces  qx,  q^,  qn  also  include 
inertia  terms.  The  following  Fourier  tranforms  of  the 
displacement  components  u,  v,  w  are  defined: 


U  4 
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J  j  u(C.n)  d£dn 


o  o 
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|  |  v(c,n) 


*m  dWn 


0  0 
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w(C.n)  ♦i?)  dCdn 
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mn 

=  cos^? 
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00) 

m/l ,  8n 

=  mr/B, 

a  =  a/R, 

b  =  b/R 

The  system  of  differential  equations  is  subjected  to  finite 
Fourier  transformations  with  kernels  (8) ,  (9)  and  (10) 
respectively. 

The  double  integrals  thus  obtained  depend  on  the 
derivatives  of  u,  v,  w.  These  integrals  must  be,  as  usual, 
integrated  by  parts  in  order  to  relate  them  to  the  transforms 
(5)*(7).  In  view  of  the  discontinuities  present  this 
standard  procedure  is  extremely  tedious.  The  method  used  in 
Part  II  of  this  report  based  on  application  of  Green-Gauss  theorem 

f 

\\  dS  -  f  Fn,  dr  ♦  f  Fn.  df  (i=l,2,3) 
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is  very  helpful  in  considerably  reducing  the  amount  of 
manipulations.  Ultimately  the  following  system  of  linear 
algebraic  equations  in  umn,  vmn»  wmn  ic  obtained,  which 
includes  integrals  involving  unknown  discontinuities  of  the 
displacements  and  of  the  slope  across  the  crack: 

!“m  *  ^  6n!  ‘  r  c“!)umn  *  T*  “m6n  vmn  '  wmn  “  Atr.,» 

‘W’n  “mn  6r>2  "  p“!)vinn  *  6r>  wmn  *  Bmn 

-  n.  “n*  '  8n  +  [k2(%.2  +  ^  *  1  ‘  IT  c“!]wmn  “  Cmr> 

m,n  =  1 ,2,. . .»  (12) 


The  quantities  A^,  Bmn,  Cmn  appearing  on  the  right  hand 
sides  represent  certain  linear  combinations  (see  Apendix  2 
for  details)  of  the  integrals  of  the  unknown  discontinuities 
of  the  displacements  Au,  Av,  Aw  and  of  the  normal  slope  A  (3w/3u) 
across  the  crack.  The  system  of  algebraic  equations  (12)  is 
solved  with  respect  to  umn,  vmn,  wmn  and  the  results  sub¬ 
stituted  into  the  inversion  formulae: 


2  .  00  00 

u(C»n)  *  —  2  u  sin6_n  +  —  2  2  u  cosa  £  sinBjf 

aF  n=l  on  n  aF  m=l  n=l  mn  m  * 

(13.1) 

?  00  ,  00  00 

v(C»h)  =  —  2  v  sina_C  +  —  2  2  v__  sina  E  cosB  y 

i6  m=l  m  ^  56  m=l  n=l  mn  m  " 

(13.2) 


w(C,n)  s  —  2  2  w  sinaE  sinB  v 

56  m-1  n*l  mn  ^  " 


(13.3) 


In  order  to  determine  the  integrals  of  the  unknown  functions 
the  boundary  conditions  at  the  crack  must  be  applied  next. 


2.  Boundary  Conditions  and  the  Characteristic  Equation 


Four  boundary  conditions  are  available  at  the  crack's 
edge  to  determine  four  unknown  discontinuities.  These  are: 

"nr,  *  0.  V».  °nn  =  0-  %  =  0  <>4> 


where  Mnn  is  the  bending  moment  in  the  direction  normal  to 
the  crack,  vn  is  the  reduced  shear  force  at  the  edge,  onn  and 
ocn  are  stress  components  normal  and  tangential  to  the  crack 
respectively.  Using  standard  expressions  for  cylindrical 
shells  (see  [8])  relating  quantities  (14)  to  the  displace¬ 
ments  u,  v,  w  and  the  formulae  for  directional  derivatives 
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one  obtains  from  (14)  the  following  conditions  at  the  edge  of 
the  crack: 


(sin2a  +  v  cos2o)  ||?-  +  (cos2a  +  v  sin2a)  |^r  +  2(1  -v)  sina  cosa 

(16.1) 

sina  [1  +  (l-v)cos2a]  |^t  -  cosa  [1  +  (l-v)sin2a]  |^r 
*  cosa  [(l-2v)s1n2a  -  (2-v)cos2a]  "  s*n  Cl-2v)cos2a 

-  (2-v)s1n2a]  L4  2  3  0 
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(16.2) 
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(I6.3) 


if  +  vtn  +  vw^  S1nct  COSa  +  (vff'  +  In’  +  sina  cosa 

•  6h  ,3u  ,  3v%  /  2  •  9  \  a 

C“  (37T  3^  (cos  a  “  sin  a)  =  0 


(16.4) 


satisfy  conditions  (16) .  In  order  to  derive  proper 
differentiation  formulae  Green-Gauss  theorem  (11)  is  applied 
again  in  the  way  shown  in  Part  I.  Consider  for  instance  the 
series  (13.1)  and  the  transform  (5).  In  order  to  find  9u/3? 

3 

Green-Gauss  theorem  is  applied  to  the  expression  u  <J> 

mn 

yielding 


f  3u  . (3) 

J  H  *mn 


♦» ds  -  -  n u 


«  ♦  (  u  ♦In’  =£  df 


where  r  is  the  "contour"  of  the  crack 

Since  <J>^  =  ct  therefore  the  above  becomes 

mn,£  m  mn 

I  f  *m  dS  *  -  %  um„  +  j  u  *m  \  dF 


The  derivative  of  (13.1)  must  be  of  the  following  form: 


a  .4-  I  Z  g  *<3> 
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where  gmn  are  Fourier  coefficients  to  be  found.  Substitution 
of  (19)  into  (18)  yields  the  desired  formula 


*  -“in  um  *  I  “  S{  dF 
f 


(20) 


Similar  approach  was  used  to  derive  expressions  for  other 
derivatives  appearing  in  eqs.  (16) . 

In  order  to  derive  conditional  equations  determining 
these  quantities,  boundary  conditions  (16)  are  applied  by 
performing  differentiation  on  eqs.  13  according  to  the  rules 
explained  in  the  previous  section.  After  lengthuy 
manipulations  this  leads  to  an  infinite  system  of  linear 
algebraic  equations  of  the  following  form 


f1  *F1U  Ui  +  F2i£  V1  +  F311  V  +  F4U  Zi  *  0 


£*0,1 


I,  <F5l£  V1  +  F6i£  Vi  +  F7i£  V  +  *  F8i£  Zi  *  0 
-i  y  w 

^  (F9i£  Vi  +  F10i£  V1  +  Flli£  Wi*  +  F12U  Zi  *  0 


£=1 . 

£=1 .... 


(F13i£  Vi  +  F14i£  V1  +  F151£  V  +  ^  F16i£  Z1  *  0 


(21) 


Coefficients  F^j£,  F2i£'*#*  ^16i£  are  composed  of  double  and 
single  series  with  complicated  general  terms,  some  involving 
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Fresnel  integrals  (due  to  the  presence  of  the  square  roots  in 
the  denominators  of  eqs.  31).  Expressions  leading  to 
evaluation  of  these  coefficiencies  are  given  in  Appendix  3 
for  the  case  when  the  tangential  inertia  terms  are  neglected. 
In  order  to  obtain  the  frequencies  of  natural  vibration  of 
the  cracked  shell  one  should  equate  to  zero  the 
characteristic  determinant  of  the  system  (26) . 

The  quantities  Aw,  Au,  and  Av  appearing  under  the 
integrals  are  not  singular  and  can  be  therefore  expanded  into 
Fourier  series.  However,  the  fact  that  the  bending  moment 
and  the  stress  are  square-root  singular  at  the  tips  of  the 
crack  makes  these  Fourier  series  very  slowly  convergent.  To 
eliminate  this  difficulty  the  procedure  adopted  in  Part  I 
will  be  generalized  here.  To  this  end  the  singular 
quantities  are  chosen  to  be  the  primary  unknows  with 
"built-in"  singularities: 


W..  sin 


iirc 

d 


/d'-c 


(22.1) 


a|  -  f  z  U,  sin  if  ♦  %  ♦  ± 

30  d  i=l  1  d  /c 


(22.2) 


A  —•  *  |  I  V.  sin  ^  +  -pz 

3n  d  i-i  1  d  v'c  &c 


(22.3) 


where  the  unknown  constants  C^,...,Cg  are  determined  using 
continuity  conditions  of  the  displacements  at  the  tips  of  the 
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crack  and  by  applying  integration  by  parts.  Variable  "c"  is 
measured  along  the  crack  starting  at  one  tip.  Extensive  but 
straightforward  manipulations  yield  finally  the  following 
expressions: 


3  00  Wi  2  ”  Wi 

4  Z  -A-  +  —  z  r*- 


1  ttv^T  i=l  1  vJd  i=1,3  1 


3  »  W.  4  “  Wi 

c,  *  —  e  r  -  -2—  z.  r5- 

^  irv^d  i=l  1  ir^  i=1 ,3  1 


'3 


!  »  (-1)\ 

ir^T  i=l  ~ 


-  a  3/2  “  Hi  ,2  HI1,  .3  r  W1, 

**  1.1  1  ‘5  21 V’  5  if,  ,3  ~] 


-  -  1 


__L  J  Hi  +  2x_  d  3/2  r  “  Wi  ,2  1  .  8  “  Wi-. 

1-1  '  M1 


!  “  (-DiVi 

ir«^  i=l  i 


2i_d  3/2  [  z  !i  (2  .  (-I)',  3  "  "“i, 

^  S=1  1  (5  i??’-  5-,fl,3^ 


00  I'W. 


00  V. 


J_  r  *i  +  3/2  r  "  Wi  ,2  1 

1=1  r  ^  d  C1=l  T  (*  ’  5? 


w. 


o  r>  ln 

T5  z  r-J 

15  1-1,3  1 

(23) 


The  only  unknown  function  which  is  not  singular  is  the  slope 
discontinuity  across  the  crack.  Tt  is  represented  by  the 
following  series: 
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f.V 
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K  ••  ' 
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The  condition  that  A  (3w/3n) be  zero  at  both  tips  leads  to  the 


following  relationships: 


-2  I  Z. 
i=2,4  1 


(25.1) 


h  ‘  ~  1  Zi 

i=3,5  1 


(25.2) 


The  unknown  functions  Aw,  A(3u/3n),  A(3V/3n)  appearing  under 


the  integral  signs  in  the  quantities  Atnn,  Bmn,  Cmn  <EcJs*  12) 
are  replaced  now  by  their  representations  (eq.  27.1*27.3, 

23) .  The  discontinuities  of  the  higher  order  derivatives  of 


the  displacements  are  eliminated  by  using  the  conditions  (14) 


referred  to  the  local  coordinate  system  Ocn  along  the  crack 


(c,n  being  directions  respectively  tangential  and  normal  to 


3  3 

the  crack) .  Thus  to  eliminate  for  instance  A  3  w/3  n)  we 


would  use  the  condition. 


A  ^  =  -  (2-v)  A  -5^- 

3(T3n 


and  then  eliminate  the  dependence  on  "c"  by  integration  by 


parts. 


Presently  all  the  quantities  appearing  on  the  right  hand 
sides  of  eqs.  12  (and,  after  inversion,  in  eqs.  13)  depend 


solely  on  four  infinite  sequences  of  Fourier  coefficients: 


U^,  V^,  ,  and  Zj;. 

The  frequencies  have  been  calculated  for  the  crack 


positioned  at  the  apex,  parallel  to  the  straight  edge  of  a 


6! 


shell  with  square  planform  for  various  d/a  ratios  and  for 
three  R/a  ratios:  R/a  =  25,  R/a  =  2.5,  and  R/a  =  1.0.  (See: 
Figs.  6,  7,  8) . 

The  first  case  is  equivalent  to  an  almost  flat  plate  and 
was  considered  for  the  purpose  of  compassion.  In  the  second 
and  the  third  cases  the  effect  of  the  curvature  is  more 
evident.  Since  no  other  data  are  available  for  these  last  c 
ases  therefore  the  frequencies  of  the  uncracked  shell  with 
identical  geometry  were  evaluated  first.  For  this  end  the 
characteristic  equation  for  a  shell  with  tangential  inertia 
effects  neglected  has  been  solved 

dmn  *  %  *  '’'frn  +  (68n  *  ^T"  ' 

+  <46n  -  2»  +  6n  -  ^  *  0  <27> 

as  well  as  the  characteristic  equation  corresponding  to  the 
case  when  all  inertial  effects  are  included.  It  was  found 
that  the  decrease  of  the  curvature  increases  the  frequencies 
and  that  neglecting  tangential  inertia  slightly  overestimates 
the  magnitudes  of  frequencies  as  shown  in  Table3. 
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Appendix  1 


The  following  is  the  list  of  symbols  appearing  in  equations  (21)  *  (24) 


X<i _ *  -2n(agcn  -  a^rn  +  a? r  ^  m  n  +  <p  m  ) 


2mn 

(3mn 


5c 


7c 


*6c 


X^__  =  n(alcn4  -  a3c  4>2m2n2  -  a2c  <fr4m4  +  a4(.n2 


.4_4 


.2_2_2 


5mn  8s 


aas  +  aQ  $  m  n  +  agc  n  -  ainil 


‘101 


X6mn  '  ♦"  (a10c*2,n2  *  a10s  "2>  -  X7mn  =  ^"(‘ios*2”2  *  a10cn2> 


.  4  4  7  7  7  4  2 

X8mn  *  2al 1 ^  m  +  2v*  m  n  +  n  - 


xlmn’  X4mn’  X7mn  are  obtai"ed  frMI  x3mn-  X2mn'  x&nn  respectively  by  replacing 
n  by  <|«n  and  a^c  by  a-s  and  vice-versa 


X9mn  *  a15s  ^  m  +  ai6s  4  m  n  +  a17s  ^  m  n  “  aigs  n  (n^_n  ) 


2_2„2 


+  ®-jqs  ♦  w  fl 


X10mn 

X13mn 

X16mn 


^•4  2  2  2  4  2 

=  4>mn(a20s  ♦  'm  +  a21$  m  n  +  a22s  n  +  a23$n  ) 

4  4  2  2  2  4 

*  4*n(®24s  ♦  m  +  a25$  ♦  n  n  +  a2gsn  + 

•  ♦mCa28cn4  +  a29c  +  2ai4c  (4»4m4-«2)] 


X1 lmn ’  X12mn*  X15mn’  X14mn  are  obtained  from  X9mn’  X10mn’  X13mn’  and 
X16mn  respectively  by  replacing  n  by  <J>m  and  aic  by  ais  and  vice-versa. 


I 


The  symbols  a-c  and  a,.^  are  listed  below  where  for  abbreviation 


s  =  sin  ip,  c  h  cos  ip 

a1$  5  -0-v2)sc,  a2s  =  (l-v)[(l-3v)c4  -  2(l-2v)c2  +  l-v]s 
a3s  -  2(1-S)(2-v)9c2(1-2c2),  a4s  =  [-1  +  (2-v)c4]s 
a5s  =  (l-v)sc4,  ags  =  -v(l-v)s3c2 
a7s  a  0”v)[-3  +  2v  +  (5-3v)sc2,  a8s  =  (l-v2)c4 
a9s  s  2(1-v)[1-(3-v)c2  +  (3-v)c4,  a]0  =  1  -2(l-v)c2  +  2(l+v)c4 

a10s  *  -2(1-v)[1-(1-v)s2]sc,  =  (l-v)s2c2 

a12  *  0-v)sc4,  a13s  =  -s4  [s4  -  vs2c2  +  (l+v)c4] 
a14s  s  -s4c(s2-3c2)  ,  a15s  =  (Hv)sZc4  (l-4c2) 
al6s  =  2£~7  +  2v  "  2(1+v)c2]s2c4 
a17s  =  "v“2(l“5v)c2  -  (7+19v)c4  +  4(l+v)c6]s2 

a18s  *  +  Z<*2  "  2vc4  +  4(Hv)c6]s2 

a19s  =  2f>(2+v)c2  +  2(1+v)c4]s4 
a20s  =  ^2(l+v)c2  -  3(5+3v)c4  +  8(l+v)c^]sc 
a2ls  *  2£“1  ”0-9v)c2  -  (7+17v)c4  +  8(Hv)c6]sc 

a22s  =  C-5  ♦  v  +  12c2  -  3(5+3v)c4  +  8(l+v)c6]sc 
a23s  =  "  4(l+v)c2  +  12(l+v)c4  -  8(l+v)c®]sc 

a24s  "  ("7+v+8c2)sc4,  a25s  =  2[-2  +  v  -  2(2+v)c2  +  8c4]sc2 


----- 


*265  =  [-1  +  v  +  2(l-2v)cZ  +  3(-3+v)c^  +8js 
a27s  =  tf-v  +  "  ^c^]s.  &28S  =  2(6-v-4cZ)sZc^ 

a29s  =  +  (7+3v)cZ  -  8c^]sZc 

The  quantities  a^  are  obtained  from  corresponding  quantities  a^s  by 
replacing  s  by  -c  and  c  by  -s. 

Also: 


u 

Tlmnk  '  2  j  <C0SV  sinV,F  sin  T  dc 

o  r 

d 

T2mnk  =  2  }  (s,nV  C0SV>F  s("  T  dc 

o  1 

d 

T3mnk  *  2  |  (sinV  sinV>F  cos  dc 


□ 

T*nnk  *  2  j  <C0SV  C0SV>f  cds  IT  d‘ 


t  _  I  ci.  3mnc  .  kirc  f  _  ( 
-  .-I  sin  j  sin  A  dc,  -  I 

0  0 


lmnk 


V  c 

4mn  kirc 

2mnk  I  S1n  d  sin -j- dc 


3mnk 


V  c 

3mn  ^  kir Cj.  i. 

COS  — g Sin  — g1  dc  f  t 


4mnk 


.....  v4mnc  knc 
COS  — g —  sin  — g—  dc 


All  these  Integrals  have  been  represented  by  closed  form  expressions.  Also 

Zlmnk*  Z2mnk  *  T~  (t3mnk  sinvlmn  +  tlmnkvcosVlmn) 
v3mn 


1 

4  (t4mnk  s(nv2mn  4  ‘zmnk  cosv2nJ 
mn 
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Where  the  upper  sign  applies  to  Z1mnk  and  where 


'lm  =  Vi  *  Vl  ■  '  Vl  -  Vl 


v3mn  *  <%  cos  *  *  6„  5in  *>d-  ■  (Vos  *  -  V*"  *>d 

Finally 

r  5/2 

Elmn-  E3mn  *  -(1/,3mn  >  ([sinvlinn  *  s,n<vl*i  *  v3nJ] 

*  Ccdsvlmn  4  cos(vlmn*  »3*,>1  s(^ 

1  <’/vUn>  C(^) 

*  [cosv^  *  cos(v2mn  *  ,%n)]  S<4^)> 


E2«n-  E*,„  ■  C/.L2)  -  2  4  *3*)] 


*  [cosvln,n  *  2  cos(vi«ntVi™n)]  S (*ZD 


lmn  3mn; 


3mn' 


4  (’/’tan)  Csi"V2»n  -  2  sin(v2«mt,*rai>^  C<^> 


4mn' 


4  [cosv2».  4  2  cos  <v  Amn  )  ^  s(^n)) 


2mn  4mn' 


4mn' 


Where  the  upper  signs  refer  to  Elmn  and  E2mn  and  where 

s(x)  |  !r^dt  c<*>  *  j  cjr 


dt 


are  Fresnel's  integrals. 


The  right  hand  sides  of  eqs.  (12),  Part  II  are: 
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(mn) 
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+  Rk2(^B2) 


A  4>3  dc  +  k2R3(l-v) 
3n  (mn) 


+  (vs2+c2)  (  Av  4>3  dc 


ix  3:  Symbols  used  in  the  series  eq.  (26) 


J 


fj  4 
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In  order  to  obtain  the  coefficients  F^,...,  F]6iJl  the  conditional 


equations  (26)  are  represented  in  the  following  form: 
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mn 
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mi 
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The  symbols  appearing  in  the  equations  (A2-1  -  A2-4)  are  listed 
below: 
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s  =  sina  c  =  cosa 
F1  =  s2+vc2 

F1  *  c2+vs2 

F2  =  1  -  2(1-v)s2c2 

Y  =  [1  +  (l-v)c2]s 

Y  =  [1  +  (l-v)s2]c 
6  =  [(2v-l )c2  +  (2-v)s2]s 

=  [(l-2v)s2  +  (v-2)c2]c 

F3  =  -s(ys2+6c2) 
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Table  1 

Evaluation  of  terms  appearing  in  c4kl,  c5kl>  cgkl 


Quanti ty 

Xlmn 

X2mn 

1 

X 

X4mn 

c 

£ 

X 

X6mn 

X7mn 

X8mn 

Substitution 

X9mn 

X12mn 

X10mn 

Xllmn 

X13mn 

X15mn 

X14mn 

X16mn 

Quanti  ty 

T3mn£ 

T4mn£ 

a5S 

a5c 

T4m0£ 

^4on£ 

allT4m0k 

allT4onk 

Substitution 

^lmn£ 

T2mn£ 

~a13c 

"a13S 

a14sT4m0k 

a14cT40nk 
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Table  2 


Comparison  of  the  Fundamental  Frequency  Factors  with  the  Results 
Obtained  in  [5]. 


-  60  - 


Table  3 


Influence  of  the  Relative  Location  (ei/a)  of  the  Crack  on  the 
Frequency  Factors  f  for  Lowest  Two  Modes  (d/a  =  0.5). 


ei/a 

f. 

f2 

Comments 

.01 

1.947 

■1 

Inaccurate 

.05 

1.924 

MEM 

.10 

1.900 

.15 

1.877 

All  Results 

o 

C«J 

• 

1.856 

4.154 

From 

.25 

1.835 

4.166 

20  X  20  Matrix 

.30 

1.818 

4.203 

.35 

1.804 

4.253 

.40 

1.794 

4.303 

.45 

1.787 

4.338 

.50 

1.785 

4.351 

Table  4 


Effect  of  curvature  and  of  tangential  inertia  on  the  two  lowest 
frequencies  of  vibration  of  an  uncracked  shell  of  square  planform 


Tanqent.  inerti 

a  not  Included 

1  Tangential  inertia  included  1 

mm 

fi 

f2 

f2 

oo 

2.000 

5.000 

2.000 

5.000 

25 

2.003 

5.000 

2.003 

5.001 

2.5 

2.269 

2.256 

1.0 

3.343 

5.112 

3.227 

5.042 
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rigure  3. 

-requency  coefficients 

Pora  square  plate  with  a  parallel  crack 
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Figure  6. 

Frequency  coefficients  fora  shell 
cracked  at  the  apex 
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Mode  2  \ 

Figure  7 

Frequency  coefficients  for  a  shell 
cracked  at  the  apex 
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Figure  8 

Frequency  coefficients  for  a  shell 
cracked  at  the  apex 
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